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Abstract. We compute the next-to-leading order QCD corrections to the "direct" part of the spin- 
dependent cross section for hadron-pair photoproduction. The calculation is performed using largely ana- 
lytical methods. We present a brief phenomenological study of our results focussing on the A"-factors and 
scale dependence of the next-to-leading order cross sections. This process is relevant for the extraction of 
the gluon polarization in present and future spin-dependent lepton-nucleon scattering experiments. 

PACS. 13.88.+e - 12.38. Bx - 13.85.Ni 



1 Introduction and Motivation 



For many years the field of QCD spin physics has been 
driven by the hugely successful experimental program of 
polarized deeply-inelastic lepton-nucleon scattering (DIS). 
One of the most prominent results has been the finding 
that quarks and anti-quarks summed over all flavors, AS, 
provide only about a quarter of the nucleon's spin, con- 
trary to naive expectations from quark models. This im- 
plies that sizable contributions to the nuclcon spin should 
come from the polarization of gluons, Ag(Q 2 ), or from or- 
bital angular momenta L q ^^ g (Q 2 ) of partons. Here, Q de- 
notes the resolution scale at which the nuclcon is probed. 

Results from fully inclusive DIS experiments are now 
supplemented by a growing amount of data from polar- 
ized proton-proton collisions at BNL-RHIC p], but also 
from less inclusive measurements in lepton-nucleon scat- 
tering [2,3,4,5]. Determining the gluon spin contribution 
is the major focus of all these experiments. The strength of 
RHIC is the possibility to study several different processes 
over a wide kinematical range which are directly sensitive 
to gluon polarization [5J: single-inclusive prompt photon, 
jet, hadron, and heavy flavor production at high transverse 
momentum px or any combination of these final-states in 
two-particle correlations. The way to access Ag(Q 2 ) in 
lepton-nucleon scattering is to select final-states which are 
predominantly produced through the photon-gluon fusion 
(PGF) process. Due to the relatively small center-of-mass 
system (c.m.s.) energy y/S available in current fixed-target 
experiments, such studies are limited to charm and single- 
or di-hadron production at moderate pr- 

To reliably determine of the amount of gluon polar- 
ization Ag(Q 2 ) entering the proton helicity sum rule, it 
is imperative to precisely map its Bjorken-x dependence 



first, in order to minimize extrapolation uncertainties in 



Ag(Q 2 ) 



Ag(x,Q 2 )dx 



(1) 



The eventual extraction of Ag(x, Q 2 ) will require consider- 
ation of all existing data through a "global QCD analysis" 
that makes simultaneous use of results for all probes, from 
pp and from IN scattering. This is the only way to effec- 
tively deconvolute the experimental information, which in 
its raw form is smeared over the fractional gluon momen- 
tum x and is taken at different scales Q. 

The basic concept that underlies the theoretical frame- 
work for high-pr processes, and any global analysis thereof, 
is the factorization theorem. It states that large-momen- 
tum transfer reactions may be factorized into long-dis- 
tance pieces that contain the desired information on the 
spin structure of the nucleon in terms of universal parton 
densities Af(x,Q 2 ), f — q,q,g, and parts that describe 
the short-distance, hard interactions of the partons. The 
latter can be evaluated within perturbative QCD. Here, at 
least next-to-leading order (NLO) accuracy is required for 
quantitative analyses to control theoretical uncertainties. 

A first such global QCD analysis of Af(x,Q 2 ) is now 
well under way [7], including all recent results from the 
RHIC experiments PHENIX and STAR [T|. However, re- 
sults on hadron-pair production from polarized lepton- 
nucleon scattering experiments [3|j4] have to be left out 
due to the complete lack of NLO computations for this 
important class of processes. Two-hadron photoproduc- 
tion is also expected to play an important role in the spin 
physics program at a future polarized lepton-proton col- 
lider, which is currently under discussion [8]. 

This paper is the first step towards a full NLO de- 
scription of hadron-pair production in longitudinally po- 
larized lepton-nucleon collisions. Here, we compute the 
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NLO QCD corrections to the "direct" part of the spin- 
dependent cross section for two-hadron photoproduction, 



l(Pl,\i)N(P N ,\ N ) - l'(P tl )H 1 (P 1 )H 2 (P 2 )X 



(2) 



i.e., where the exchanged photon is at low virtuality and 
interacts as an elementary particle with one of the partons 
of the nucleon N. The Pi in © are the four-momenta 
of the observed leptons and hadrons, X contains all the 
additional hadronic activity not observed in experiment, 
and the A<j denote the helicities of the interacting lepton I 
and nucleon N. 

Of course, an immediate complication arises here, as 
the direct part on its own is no longer a well-defined quan- 
tity beyond the leading order (LO) approximation. This is 
due to kinematical configurations with a collinear splitting 
of the photon into a qq pair which need to be factorized 
into the photon structure functions appearing in the "re- 
solved" part of the cross section. This is well-known [9] 
and expresses the freedom in the factorization procedure 
such that only the sum of direct and resolved contribu- 
tions is independent of theoretical conventions. Neverthe- 
less, we will concentrate in this work on the direct part of 
the polarized two-hadron production cross section which 
is technically already rather involved. 

This is because we perform the calculation using large- 
ly analytical methods. Calculations of this kind were pi- 
oneered in the unpolarized case for photon- hadron [15] . 
photon-photon [TT], and photon-charm [TJ] correlations 
quite some time ago. In the polarized case only a calcu- 
lation for double-photon production [13] exists so far. To 
keep the computations tractable, we chose to present the 
results in terms of the transverse momentum and rapid- 
ity of "hadron one" (Hi), Pt,i and yi, respectively, and a 
variable Q3HT] 



ZH 



TA 



T,2 



p2 



(3) 



which contains some information about the kinematics of 
"hadron two", H 2 , but not including its rapidity y 2 . Nu- 
merical evaluations of the triple differential cross section 
are thus limited in that experimental cuts on the rapid- 
ity of the second hadron H 2 cannot be implemented. The 
introduction of zh in the analytical calculation is essen- 
tial to keep certain singular configurations at bay [101111] , 
for instance, the case when the two hadrons are produced 
collinearly, as will be discussed in more detail below. 

Despite the fact that our calculation is not complete in 
the sense discussed above, we strongly believe our results 
to be very important, both theoretically and phenomeno- 
logically. On the one hand, our results will serve as a check 
on more versatile calculations in the future using com- 
bined analytical and Monte Carlo techniques which we 
are pursuing at the moment [14] along similar lines as in 
the unpolarized case [15] . These studies [14] will include 
also the spin-dependent resolved photon contributions at 
NLO accuracy. On the other hand, it was demonstrated 
in recent a LO study [16] that the direct photon part is 
responsible for the main features of the experimentally rel- 
evant spin asymmetry and its sensitivity to the polarized 



gluon density at fixed-target experiments like COMPASS 
and HERMES. The resolved photon part is non- negligible 
though, but merely leads to a roughly constant shift of 
the spin asymmetries. We also believe that our numerical 
studies of the relevance of the NLO corrections and theo- 
retical uncertainties due to variations of the factorization 
and renormalization scales already give a good indication 
of what to expect from a full NLO calculation in the fu- 
ture. 

The paper is organized as follows: in Sec. 2 we present 
the details of the calculation of the NLO QCD corrections 
to the direct part of spin-dependent two-hadron photopro- 
duction. Section 3 is devoted to a brief numerical evalua- 
tion of our results, focussing on the relevance of the NLO 
corrections and the residual scale uncertainties at fixed- 
target kinematics. Section 4 contains the conclusions. In 
the Appendix we collect some additional details of the 
calculation. 



2 Details of the Calculation 

2.1 General Framework 

The process we want to consider in the following is the in- 
clusive production of a pair of hadrons Hi H 2 in collisions 
of longitudinally polarized leptons and nucleons with four- 
momenta as specified in Eq. ([2]) . Both hadrons are required 
to be at high transverse momentum. As mentioned above, 
we consider only the direct part of the cross section, where 
the exchanged photon interacts as an elementary particle. 

Since we want to perform the NLO calculation using 
largely analytical methods, we are limited to observing a 
hadron Hi with transverse momentum Pt,i and rapidity 
j/i, together with hadron H 2 in the opposite hemisphere, 
its transverse momentum vector Pt.2 constrained by zh 
defined in Eq. ([3]), but otherwise unspecified kinemat- 
ics. Assuming, as usual, factorization, we may then write 
the NLO expression for the corresponding spin-dependent 
cross section as a convolution of the non-perturbative par- 
ton distribution and fragmentation functions and the hard- 
scattering of the partons 



dAo HlIi2 
dPr,idyidzH 



aa ++ 



do 



H1H2 



dPr,idyidzH dPr,idyidzH 



2P. 



TA 



S 



E 



dzi 



k Ji-v+vw z i 



dv 



(4) 



dw 



dz 



—Af! y (x l ,ii f )Aff(x N ,f, f ) 



v(l — v) J vw_ w 



dv 



5(1 - w)S(l -z) + 



Os(pr 

2tt 



dAa (1) -,y 

X — : (s,V,W,fJ,f,fl f ,fJ, r ,z) 



dvdwdz 



xD^(zi,^)D^(z2,^' f ) 



(5) 
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The subscripts "++" and "H — " in J3J denote the settings 
of the helicities of the incoming lepton and nucleon. We 
have introduced the standard hadronic invariants 

S=(P l + P N ) 2 , T=(P l - Pi) 2 , U=(P N - P) 2 , (6) 



V = 1 + - w = — — 



(7) 



with four- moment a specified in and their partonic 
counterparts 

s = (p 7 +Pif, t = (p 1 -pj) 2 , u = (pi -Pj) 2 , (8) 



v = 1 + -, w 



s+t 



(9) 



Neglecting the masses of all particles one finds the follow- 
ing relations among the variables in Eqs. ©-© 



a = x lXN S, t = fT,u= z*-U, 

_ VW _ l-V 

I vwz± ' ^ [T— vjzl 



(10) 

(11) 



with x e [xn] the fraction of the longitudinal momentum of 
the lepton [nucleon] taken by the quasi-real photon [parton 
i]. In addition, V and W in (J7J) are determined by the 
observed hadron H\ and the lepton-nucleon c.m.s. energy 
squared S: 



i/ i Pt ^ - 
V = 1 =-e 



l , w = 



p2 



SV(l-V) 



(12) 



In ((5l) we have also introduced the partonic counterpart 
of zh, defined via 



PT,j ■ PT.k 

Ph 



£l 

Z2 



(13) 



with px.j and pt,h the transverse momenta of the final- 
state partons i and j producing hadrons Pi and P 2 , re- 
spectively. z 12 are the momentum share that the hadrons 
Pi. 2 inherit from its parent partons j, k in the hadroniza- 
tion process. The latter is modeled by non-perturbative 
functions D^' 2 (21,2, /-*/) describing the collinear fragmen- 
tation of the partons j and k into the observed hadrons 
Pi and P2, respectively. 

The Af N (xN, pf) in (O are the spin-dependent parton 
distribution functions, defined as usual by 

Atfixs^f) = f t N {+) (x N ,H) - f«-)(*N,lif) ■ (14) 



The subscript + [— ] in Eq. ([14]) indicates that the parton's 
spin is aligned [anti-aligned] with the spin of the parent 
nucleon N. 

As we only focus on the direct photon case in our 
calculation, AfUxi,fif) in ([5]) coincides with the spin- 
dependent Weizsacker- Williams equivalent photon spec- 
trum, which reads [17] 



2tt 



1-(1 



In- 



xi 



+ 2mfx 2 ( 



) 2 

v max 



1 - Xl . 

m 2 x 2 ■ 



(15) 



with mi the mass of the lepton. It describes the collinear 
emission of a photon with low virtuality Q, less than some 
upper limit Q max determined by the experimental condi- 
tions. The non-logarithmic pieces in (| 15(1 result in a small 
but non-negligible contribution in case of muons. 

The sum in |JSJ runs over all possible partonic channels 
7? — > jkX, with dAa^ , and dAa ■ , v the associ- 
ated LO and NLO longitudinally polarized partonic hard- 
scattering cross sections, respectively. They are defined in 
complete analogy to Eq. ((4]) and have been stripped of 
trivial factors involving the electromagnetic coupling a em 
and the strong coupling a s (p r ) evaluated at renormaliza- 
tion scale fi r . As indicated in ([5]), starting from the NLO 
level, the subprocess cross sections will explicitly depend 
on n r , as well as on the scales fif and jjb'f of the par- 
ton distribution and fragmentation functions owing to the 
factorization of initial and final-state collinear singular- 
ities to be discussed below. The calculation of the NLO 
dAa^^ kx / dvdwdz is the main purpose of the remainder 
of this paper. 

A computation with largely analytical methods be- 
comes feasible [T0 l [TT ]rr2iri3] thanks to the introduction of 
the variable z (or Zjj) describing P 2 . The price to pay is a 
limited control of the kinematics of hadron P 2 , most no- 
tably its rapidity y 2 - The main virtue of z is that when in- 
tegrating over phase-space certain singular configurations 
of two partons can be easily avoided. For instance, mo- 
menta Pj parallel to pu corresponds to negative values 
of z. In this case, the factorized expression (|4]) for the 
cross section is incomplete and contains uncanceled poles 
which would require the introduction of additional non- 
perturbative functions describing the simultaneous frag- 
mentation of a single parton into two hadrons. Situations 
where hadron P 2 is parallel to the direction of the incom- 
ing photon are characterized by z = 0. 

For the phase-space integrations of the matrix ele- 
ments we therefore restrict ourselves to 



z > z n 



>0, 



(16) 



equivalent to the condition that the two hadrons are pro- 
duced by partons in opposite hemispheres. The final ex- 
pressions for the subprocess cross sections will contain 
mathematical distributions in 2, i.e., 5(1 — z), 1/(1 — 
etc., in addition to similar functions in the variable w, 
which are already present in the NLO computation of 
single-inclusive hadron spectra [T81IT9] . The analytical re- 
sults will therefore be rather involved and lengthy. Some 
technical details may be found in the Appendix. 

Finally, we note that corresponding expressions for 
spin-averaged cross sections are straightforwardly obtained 
by replacing all polarized quantities in this subsection by 
their unpolarized counterparts. To the best of our knowl- 
edge, NLO corrections have not been computed with an- 
alytical methods in the unpolarized case either. We will 

therefore provide results also for da^_^ - kx /dvdwdz. 
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2.2 LO Contributions 

As we consider only the direct part of the photoproduc- 
tion of two hadrons, there are only two partonic chan- 
nels contributing to ([5]) in the lowest order approximation: 
photon-gluon fusion, 



19 -* , 

and the QCD Compton process, 
79 -> 19 ■ 



(17) 



(18) 



Since at LO s + t + u = and pxj — ~Pr,k, the partonic 
cross sections are <5-f unctions in w and z, i.e., 



dvdwdz 



dAa {0) 



dv 



5(1 - w)5(l - z) (19) 



as indicated in (O. The PGF process is symmetric under 
exchange of q and q, and the result for jq — > can be 
obtained by replacing w by 1 — v. Explicit expressions for 
dAa^^^ can be found, e.g., in [18] . Phenomenological 
studies based on LO results have been performed, for in- 
stance, in [20] and, most recently, in [16] . The unpolarized 
counterparts da^^^ k are given, e.g., in [21] . 



2.3 Computation of NLO Corrections 

At NLO, three different types of contributions have to be 
considered and evaluated: 

1. the interference of the tree- level amplitudes for the pro- 
cesses (fT7| and (|18[) and the virtual, one-loop correc- 
tions to them; 

2. the real gluon emission corrections to the tree-level 
processes, i.e., jq — > qgg and — > qqg; 

3. genuine NLO processes, i.e., ^q — * q'q'q, jq — > ggq. 

To account for singularities one encounters when cal- 
culating the loop diagrams or when performing the phase- 
space integrations for the unobserved parton, we use di- 
mensional regularization, where space-time is extended to 
n = A — 2s dimensions. To project onto definite helic- 
ity states for the incoming parton and photon, we adopt 
the standard HVBM prescription [25] to define 75 and 
the Levi-Civita tensor in n dimensions. The relevant n- 
dimensional partonic hard-scattering matrix elements be- 
fore integration over phase-space are the same as for single- 
inclusive particle or jet production and hence well known. 

The virtual corrections can be found, e.g., in [18j and 
[21] in the polarized and unpolarized case, respectively. 
Since they resemble the two-body final-state of the LO 
result in (TT9"]) , their contribution to the 0(a em a?) correc- 
tions is proportional to 6(1— w)S(l — z). The polarized and 
unpolarized NLO matrix elements in n dimensions with a 
thrce-parton final-state can be taken from Refs. [ 1 8 , 2 1 J as 
well. Integrating them analytically over the phase-space 
of the unobserved parton is, however, much more involved 
than for single-inclusive hadron production and special 
care has to be taken. 



2.4 Phase-Space Integration 

We consider a generic "2 — * 3" photoproduction process 

l(Pj)i(Pi) -> j(Pj)KPk)l(pi) (20) 

contributing to d(A)a^^- kx in Eq. ([5]). Partons j and k 
shall produce the two observed hadrons in the fragmenta- 
tion process. Parton I remains unobserved and hence has 
to be integrated over the entire phase-space. The phase- 
space of one of the observed partons, say fc, is constrained 
by z but otherwise integrated. We shall perform all inte- 
grations analytically and express the result in terms of the 
variables v, w, and z. 

Starting from the definition of the three-particle phase- 
space in n dimensions 

dP S 3 = [ f p * ^\ Pk J" 131 S(p?)6(pl)6(pf) 

J (27T)"- 1 (27T)"- 1 (27T)"- 1 K ] ' K k> UU 

(27T) n S^( Pl + Pl ~ Pj - Pk ^ Pl ) (21) 

one proceeds at first along the same steps as for a one- 
particle inclusive final-state [23] , arriving at the well-known 
result 

2e 



dPS, 



Air 
s 



,l-2e 



(1 



l dv 



(47r) 4 r(l - 2e) 
x[w(l-w)]- £ dw J d6 1 d6 2 (sm8 1 ) 1 - 2£ (sm6 2 y 2£ .(22) 

To proceed, we parametrize the momenta in (|20p in the 
c.m.s. frame of partons k and I, 



(l,p x ,sin0i cos 2 , cos 9i, Pk ) 
( 1 , -p x , - sin 0i cos 2 , - cos 0i , - Pk ) , (23) 



P k 



Pi = 



2 



where s k i = (pk +Pi) 2 = sv(l — w). pk denotes the (n — 4)- 
dimcnsional components of Pk and Px is arbitrary. The 
other three momenta can be chosen to have non- vanishing 
spatial components only in the y- and z-directions, and 
explicit parametrizations can be found in the Appendix. 
The variable z in (|13p is introduced as 

z = m- Pk (24) 

by defining an auxiliary space-like vector jTOlfTT) 

_ P ~ f u+ Pi t +pjS 



tu 

and using the identity 

1 = J dz5(z — m ■ pk) 

to perform the integration over Q\ in 
arrives at 

2e 



(25) 
(26) 

One finally 



dPS* = 



(47r) 4 r(l - 2e) 



47T 

s 



: (l-v)- £ dv 



x[w(l — w)] £ dwdz2 



w(l — v) 



1 — vw 
1 — w + 4u>(l — v)z(l — z) 



1 



vw 



d9 2 sin" 25 2 , (27) 
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in agreement with the result given in [lOpilj. Further in- 
tegration over 02 depends on the various combinations 
of scalar products of parton's momenta appearing in the 
hard-scattering matrix elements. Structures with a com- 
plicated dependence on 2 have to be decomposed into a 
basic set of calculable integrals using momentum conser- 
vation and extensive partial fractioning. An extensive list 
of basic integrals can be found in [TUHmfT21fT5] and need 
not be repeated here. We only note that after integration 
over 62 one ends up with "plus-distributions" both in w 
and in z. We refer to the Appendix and Refs. [TUIITTlfTSl 
113] for a more detailed discussion. 

Equation (|27|) is sufficient for integrating the unpo- 
larized matrix elements and most of the terms in the 
longitudinally polarized ones. In the latter case, an ex- 
tra complication arises, however, from contributions pro- 
portional to the (n — 4)-dimensional components pk in 
(|23|) . the so-called "hat-momenta". We encounter terms 
in the polarized matrix elements which depending on p\ 
and require modifications to (|27p . since in its derivation 
we have assumed that the (n — 4)-dimensional part can 
be trivially integrated. Upon a careful re-examination of 
the steps leading to (|27| , we find that these contributions 
lead to a modified phase-space formula given by 



dPS 3 = (-£)- 



(4^)4^(2 - 2e) 
x (1 — w)[w(l — w)]~ e dw dz 
1 — w + 4w(l — v)z(l — z) 



47T 

s 



2= 



,2-2e 



(1 -v)- e dv 



w(l — v) 
1 — vw 

1-e 



J d0 2 sin 2 - 2e 02 , (28) 



1 — vw 



in agreement with the result given in Ref. 13J3- As is 
explicit in (|28p , all contributions stemming from p 2 , are of 
0(e) as they should. Nevertheless, they can lead to finite 
contributions in the limit e — > whenever they pick up a 
1/e pole in the remaining phase-space integrations. 



2.5 Factorization and Final Results 

Adding the renormalized virtual corrections and the real 
contributions, all infrared singularities cancel out, includ- 
ing all 1/e 2 terms. The remaining 1/e singularities are of 
collinear origin and arise when the momentum pi of the 
unobserved parton I becomes parallel to any of the other 
parton momenta. Singular configurations related to the 
initial-state are absorbed at a factorization scale [if into 
the definition of the parton densities. Similarly, final-state 
mass singularities are factorized at a scale [i'f into the 
bare fragmentation functions. This is the essence of the 
factorization theorem. A special role play the singulari- 
ties from a collinear splitting 7 — > qq, which are absorbed 
into the photon structure functions. Due to the freedom in 
choosing the amount of finite pieces subtracted along with 



We note that some of the equations in the Appendix of 
Ref. [13] contain obvious misprints. 



the pole terms, only the sum of direct and resolved pho- 
toproduction cross sections is independent of theoretical 
conventions at NLO and beyond. We choose the common 
MS scheme throughout. 

As already mentioned in Sec. 2.1, the main virtue of in- 
troducing the variable z and demanding z > is to avoid 
certain singular contributions which are beyond the factor- 
ized framework outline here: (a), when the two produced 
hadrons are collinear, corresponding to negative values of 
z and (b), when hadron H2 is produced parallel to the 
direction of the incoming photon, characterized by z = 0. 

The factorization procedure is performed in the stan- 
dard way [23] by adding an appropriate "counter cross 
section" dAa^ act to each partonic subprocess. At NLO, 
with partons j and k being observed as hadrons Hi and 
H2, there are in principle four possible collinear configu- 
rations, and dA<jf act schematically reads 



1 dAa fact 



sv dvdwdz 



"2tt 



-AH mj [w,/j,f]- 



dAa 



mi — >jk 



[ws, v]S(l — z) 



1 



s(l — vw) 
dAo. 



AH r , 



7m — yjk 



dv 

1 



1-v 
1 — vw ' 
1 - v 



s(l — v + vw) 
dAa F 



■yi — >mk 



dv 



5(1 -z) 

1 — vw 

H jm [1-V + VW, [J,' f ] 

S(zi - z) 



1 — V + VW 



1 H W n n dAa r^jm 

tlkm[Z,[lf\- 



dv 



s,v 



x0(l -z)S(l-w) 



(29) 



The dA<7 e ab ^ cd \(^s , £] / 'dv are the n-dimensional 2^2 cross 
sections for the process ab — > cd to be found in the HVBM 
scheme in [24] . These cross sections are evaluated at some 
shifted kinematics denoted by {(s, £], since the collinear 
parton j takes away a certain fraction of the available 
momentum. Furthermore, 



(A)H ab ( K ,fj,) = ~(A)P ab (K)(^ 



+ (A)hab(K), (30) 



where 1/e = 1/e — + m47r in the MS scheme and Z\ = 
1/(1 — v + vw). In (|30|) the (A)P a b(z) denote the usual un- 
polarized (polarized) one-loop splitting functions in four 
dimensions. Note that the unpolarized H a t contributes to 
the factorization of final-state singularities since we do 
not consider the production of polarized hadrons. The 
functions (A)h a b(K) represent the freedom in choosing a 
factorization prescription, and they all vanish in the MS 
scheme, except for Ah qq (n) = —16(1 — k)/3 [5S]. Need- 
less to say that an equation similar to (|3"0|) holds in the 
unpolarized case. 
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For all subprocesses, the final polarized (and unpolar- 
ized) partonic cross sections at NLO accuracy in ([5]) can 
be schematically cast into the following form 



dvdwdz 



AK x {v, w)S(l - z) + AK 2 (v, w)5(z - z%) 

, s Qi\-z) „ , ,6(zi-z) 
(1 - Z)+ [Zi - z)+ 

ATS , ^Q( Z - 1 ) ATS , \ ( Z - Z l) 

(z-l)+ {z-zi)+ 

ATS , \ fH 1 ~ Z ) 

AK 7 (v,w) ( K 1 _ z ! 
AK 8 (v,w,z). (31) 



The coefficients AKi, i = 1,...,8, contain, in general, 
distributions in uu, and they can be decomposed further 
afl 



AKi(v,w) = Aki{v)5(l -w) + Ak 2 (v) 



1 



(l-w) + 



+ Ak 3 (v) J } ) + Aki(v, w) . (32) 

For brevity we have suppressed any dependence on the 
renormalization and factorization scales in (I31|) and f|32[) . 
Again, expressions similar to (|3T|) and (I32|l hold for each 
unpolarized NLO partonic subprocess. 



3 Phenomenological Results 

We now turn to a brief numerical study of our results 
focussing on the relevance of the NLO corrections and 
the residual scale uncertainties for both the polarized and 
unpolarized cross section. We postpone a detailed phe- 
nomenological study to a future publication [14] until the 
resolved photon contribution becomes available as well. 
For our studies here, we choose the kinematical setup 
of the COMPASS experiment, which scatters a beam of 
polarized muons with an energy of — 160 GeV off 
deuteron in a polarized 6 LiD solid-state target, correspond- 
ing to a lepton-nucleon c.m.s. energy of VS ~ 18 GeV. 

The results we show will be differential in the trans- 
verse momentum Pt,i of hadron Hi and integrated over 
the angular acceptance of the COMPASS experiment, i.e., 
covering scattering angles of less than 180mrad. Using 
y = — lntan(6*/2) this straightforwardly translates into 
a lower bound on the pseudo-rapidity y± for hadron Hi. 
Kinematics dictates the upper bound, depending on the 
hadron's transverse momentum Pt,i- Recall that we can 
not control the rapidity of hadron H 2 in our analytical cal- 
culation, which in turn implies that it may end up outside 
the acceptance of COMPASS. The range of the transverse 
momentum vector Pt,2 of H 2 is restricted by demanding 



zh > 0.4, with zh defined in Eq. ([3]). The momentum dis- 
tribution of the quasi-real photons radiated off the muons 
is described by the Weizsacker- Williams spectrum given 
in (Ti"5|) . with mi = m M and Q, 2 nax = 0.5 GeV 2 . The pho- 
ton's momentum fraction xi is restricted to be in the range 
0.1 < xi < 0.9. 

In the computation of the LO and NLO unpolarized 
cross section we use the LO and NLO CTEQ6 parton den- 
sities |26j and strong coupling a s , respectively. In the po- 
larized case, we use a special set of the GRSV helicity- 
dependent parton densities |27j , characterized by a small 
negative total gluon polarization of Ag = —0.15 at the low 
input scale of GRSV. A small gluon polarization, either 
positive or negative, is indicated by all presently available 
data sensitive to Ag(x, [if) [1,2,3,4], at least in the range 
of momentum fractions predominantly probed by these ex- 
periments, which roughly amounts to 0.05 < x < 0.2. We 
note that for a set with a small positive gluon polarization, 
like the "standard scenario" of GRSV [57] , which is also in 
agreement with current data, one encounters strong can- 
cellations between the contributions from PGF and the 
Compton process, leading to sign changes in the polarized 
cross section. This makes it rather awkward to display 
the results for the NLO corrections and the scale depen- 
dence we are interested in here. Hence, for our purposes 
we resort to the choice of Ag = —0.15. In the forthcoming 
publication [14] , we will discuss in detail the sensitivity of 
two- hadron photoproduction to Ag{x, [if). 

To model the hadronization of partons j and k into 
the observed hadrons Hi and H 2 , we use the novel set 
of fragmentation functions of DSS [28J. This new set is 
based on a first global QCD analysis of inclusive hadron 
spectra in electron-positron annihilation, DIS multiplici- 
ties, and hadron-hadron scattering and known to describe 
hadronization fairly well also in the energy range relevant 
for COMPASS [28]. Since COMPASS does not identify 
different hadron species [1] and measures only the sum 
of charged hadrons, we use the appropriate LO and NLO 
sets of DSS [28] for all our calculations. 

Figure[T]shows our results for the Pt.i -differential cross 
section for the polarized and unpolarized photoproduction 
of a pair of charged hadrons at LO and NLO accuracy at 
COMPASS. We have set all renormalization and factoriza- 
tion scales in ([5]) equal to twice the transverse momentum 
of hadron Hi . The sum of the transverse momenta of both 
hadrons might be a better motivated choice, but we have 
no control over Pt j2 within our analytical calculation. The 
so-called "if-factor" , defined as the ratio of NLO to LO 
unpolarized (polarized) cross sections 



K 



d(A)a 



NLO 



d(A)<7 



LO 



(33) 



2 The coefficients (A)Ki are too lengthy to be given here, 
but are available upon request from the authors. 



is depicted in the lower panel of Fig. [T] The computed 
QCD corrections are such that the NLO results are be- 
low the LO estimates in the entire range of Pr,i shown 
in Fig. [1] They appear to be more sizable in case of the 
polarized cross section. The observed difference of the un- 
polarized and polarized fC-factors clearly indicates that 
NLO corrections are relevant also for studies of double- 



C. Hendlmeier et al.: Next-to-leading order QCD corrections to spin-dependent hadron-pair photoproduction 



7 



d(A)a h+ + h 7dP T1 [pb/GeV] 




polarized 



j a NLO,,. LO 

dAa /dAa 



1.5 



2.5 

P T1 [GeV] 



Fig. 1. Upper panel: unpolarized and polarized photoproduction 
cross section for a pair of charged hadrons, fid — > (h + + h~)(h + + 
h~)X, at LO (dashed line) and NLO (solid line) accuracy using 
COMPASS kinematics. The dotted curve labeled "LO" refers to 
a LO calculation using NLO parton densities and fragmentation 
functions (see text). The lower panel shows the corresponding 
ratios of NLO to LO cross sections (K-factor). 



10 
10' 
10' 
10 



10* 



10 



10 



do h++h 7dP T j [pb/GeV] 



NLO 




2.5 

P T1 [GeV] 

Fig. 2. Scale dependence of the LO and NLO unpolarized (a) 
and polarized (b) cross sections for fid — > (h + +h~)(h + + h~)X 
shown in Fig. [TJ All scales are varied simultaneously in the range 
v2Pt,i < fir = fJ-f — M/ 5= 2\/2Pt.i- Solid lines correspond to 
the choice where all scales are set to 2Pt,i- All LO computations 
have been rescaled by a factor 0.01 to better distinguish them 
from the NLO results. 



spin asymmetries, All = dAa /da, as they do not cancel 
in the ratio. The contrary is often assumed in analyses of 
spin asymmetries. 

We wish to make two further remarks about the re- 
sults shown in Fig. [TJ Firstly, finding if-factors smaller 
than one is not a result of the NLO corrections to the 
hard-scattering partonic cross sections. It mostly stems 
from the difference between the LO and NLO parton dis- 
tribution and fragmentation functions, in particular the 
latter. The LO and NLO sets of DSS show pronounced 
differences, mainly because the LO fragmentation func- 
tions try to make up for the often large NLO corrections 
in some of the fitted cross sections, see [35] for details, 
which are missing in a consistent LO analysis. The effect 
of the fragmentation functions is illustrated in Fig. [TJ in 
the unpolarized case by the curves labeled "LO". They 
refer to a calculation using LO matrix elements, but NLO 
parton densities and fragmentation functions. This clearly 
demonstrates the inadequacy of LO results. At best, they 
can serve as a rough estimate, but they are insufficient for 
any quantitative analysis. Very similar observations can 
be made in the polarized case. Here, a "LO"-type calcula- 
tion leads to A"- factors ranging from 1.1 to 1.8 (not shown 
in Fig. [TJ for clarity). Secondly, we note that the details 
and size of the NLO corrections in the polarized case de- 



pend significantly on the still largely unknown gluon po- 
larization. As already mentioned, there can be strong can- 
cellations between different subprocesses, leading to sign 
changes in the polarized cross section. In their vicinity 
large NLO corrections are in general inevitable. 

As an estimate for the sensitivity of the computed cross 
sections to the actual choice of scales fj,f, p!p and fj, r in 

(J5]), we vary them simultaneously in the range v2Pt,i < 
fj, r = fXf = fj!f < 2-\/2Pt,i- We note that in principle all 
scales can be varied independently. The shaded bands in 
Figs. [2] (a) and (b) indicate the resulting residual scale 
uncertainty of the unpolarized and polarized photopro- 
duction cross sections, respectively. We find that the NLO 
results show somewhat reduced theoretical ambiguities, in 
particular, for the polarized cross section. However, simi- 
lar remarks as above apply also here. Scale ambiguities in 
the polarized case depend on the details of the helicity- 
dependent parton densities and on possible cancellations 
among different subprocesses. One also has to keep in 
mind that the results so far only include the direct pho- 
ton contribution to the photoproduction cross section. It 
remains to be seen what the effects of the resolved photon 
contribution are. We will address all these questions in a 
forthcoming publication. 
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4 Summary and Conclusions 

To summarize, we have computed, with largely analyti- 
cal methods, the NLO QCD corrections to the direct part 
of the spin-dependent cross section for hadron-pair pho- 
toproduction. This is the first step towards a full NLO 
description of this process, which plays an important phe- 
nomenological role in the determination of the gluon po- 
larization in polarized lepton-nucleon collisions studied by 
HERMES and COMPASS at present and, hopefully, at 
higher c.m.s. energies at some facility like eRHIC in the 
future. 

We find that the NLO corrections are essential for any 
study of double-spin asymmetries. They are sizable and do 
not cancel in the ratio. Theoretical ambiguities due to the 
choice of the arbitrary renormalization and factorization 
scales are somewhat reduced if NLO corrections are taken 
into account. 
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Appendix 

Parametrization of Momenta 

To perform the phase-space integration over the angle #2 
for all 2 — > 3 subprocesses ji — > jkl analytically, we choose 
to work in the c.m.s. frame of the observed parton k and 
the unobserved parton I. Their momenta are parametrized 
in Eq. ([23]) . All partons are assumed to be massless. The 
remaining three momenta are chosen in such a way that 
they have non-vanishing components only in two spatial 
directions: 

sv 

P 7 = n , — (l,0,sm^,cosVv-.)i 
s(l — tiro) 

Vi = „ , — (l,0,-smV,cosW--)) 

2,/Sfci 

Pi = S(1 :!±^ ) (l, 0, sin <//, cos <&',...). (34) 



The ellipsis in ([3"4"| denote zeros in (n — 4)-dimensional 
components, Ski = sv(l — u>), and where 



cost/' 
simp 
cos ip' 
sin ip 1 



w(l — v) 



1 — w 
1 — vw 
1 + v — vw 
1 — V + vw 
1 — v — vw 



1-v 



vw 



COS')/' , 

sin^ , 



(35) 



with v and w defined in 

The five momenta can be used to define ten different 
scalar products or Mandelstam variables. Due to momen- 
tum conservation, p 1 +pi = pj +pk+Pi, only five of the ten 
scalar products are independent. We make extensive use 
of all the relations among different Mandelstam variables 
to reduce the NLO matrix elements to a form amenable 
to analytic integration. 

In the parametrization (|23|) and (|34|) . the auxiliary vec- 
tor to, introduced in (|25|) . reads 




,0,0, 



f 



1 — w 



(36) 



Plus-Distributions 

In the analytic calculation for two-hadron production one 
encounters not only the usual plus-distributions in the 
variable w [23] . but also a host of distributions in z, the 
partonic counterpart of zh, as indicated in Eqs. ([29)) and 
(f3"Tj) . For completeness, we collect here only some of the 
definitions and identities, more details can be found in 

Refs. [T0im[l2[13]. 

With the common definition of plus-distributions via 
some arbitrary test function f(z) [23] j one nas 



( z i - z )+ Jo Zl-Z 



dz- 



(z- Zl)- 



dz f{z) - f[zi \ (37) 



z — z\ 



where z\ = 1/(1 — v + vw) and z max the upper kinematical 
limit for the z-integration in ([5]) , 



w(l 



1 — vw 



-1/2 



(38) 



Analogously, the definition of other plus-distributions is 
obtained by replacing z\ by 1 in (|37j) . 

Note that in the case w = 1, the upper integration 
limit (|38[) of z becomes z max = 1, and the distributions 
at z = zi coincide with the distributions at z = 1. Also, 
in the region z > 1, the variable w is limited to ro < 1, 
and distributions 1/(1 — w) + reduce to ordinary functions. 
Finally, the lower limit z m ; n for the integration over z in 
([5]) is a function of zh and the other integration variables, 
giving rise to additional logarithmic contributions, e.g., 



1 



1 



+ 5(zi - z) ln(zi - z min ) , (39) 



(z\ - z) + (zi ~ z) Zmla 
where the new distribution is defined by 



1 - m. d*=r f ^- f ^ dz . (40) 

lin (zi - z) Zmin J Ziain zi- z 
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